Solitons in interacting Dicke models of coupled cavities with two-level systems 
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We consider an array of coupled optical cavities, each containing a multi-atom ensemble. We show that the 
nonlinearity inherent in the cooperative dynamics of the atoms in each ensemble coupled to the respective cavity 
field allows for the formation of solitary waves. Such a prediction can be tested in state-of-the-art semiconduct- 
ing photonic-crystal microcavities with embedded impurities. 
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Very recenty, a considerable amount of interest has been di- 
rected towards registers of coupled atom-cavity systems as po- 
tential candidates for the observation of cooperative phenom- 
ena in strongly correlated many -body systems [1]. Superfluid- 
to-Mott-insulator quantum phase transitions and glassy polari- 
tonic phases have been theoretically predicted in arrays of mu- 
tually coupled qubit-resonator devices. At the basis of these 
interesting effects is the nonlinear character of the so-called 
polariton, a combined state of an atom and a photon. More- 
over, the flexibility of these systems, which makes them ex- 
ploitable quantum simulators, can be used in order to realize 
effective multiple-spin dynamics useful for the purposes of 
quantum information processing |2fl. 

However, the possibilities offered by arrays of cavities in- 
teracting with atoms are enormous and intriguing. For in- 
stance, photonic Mott insulators can be achieved from atom- 
mediated nonlinear photon-photon interactions y[]. In gen- 
eral, the combination of strong intracavity atom-photon in- 
teractions and intercavity tunneling is greatly suitable for 
nonlinear-optics purposes. 

Here we further explore these possibilities and report for the 
first time the formation of solitary waves (or solitons) J3, IH 
in a periodic array of Dicke systems It is believed that 
solitons, which are localized nonlinear waves characterized 
by striking stability properties with respect to dynamical per- 
turbations (such as collisions with other solitons) represent, 
together with squeezing, one of the most remarkable manifes- 
tations in nonlinear optics |7]. The situation we consider is 
different from a classical array of media with Kerr nonlinear- 
ity as, technically, we do not expand the atomic polarization 
as a power-series of the field or assume fast response of the 
atomic system relative to the field. Differently, we treat the 
dynamics of both cavity fields and atomic ensembles on equal 
footing, our source of nonlinearity being intrinsic to each en- 
semble of atoms. The identification of soli tonic behavior in a 
system of current theoretical as well as experimental interest is 
an important step forward in the grounding of coupled-cavity 
registers as reliable and exploitable quantum simulators. We 
show that, by embedding a small number of two-level systems 
within a cavity constituting one site of an array, it is possible 
to achieve a degree of nonlinearity which, combined with the 
periodicity of the cavity arrangement, allows for the formation 
of solitons. Our results are derived by means of a simple tech- 



nique permitting the consideration of any desired/appropriate 
order of nonlinearity. Furthermore, although for simplicity we 
consider a linear array of coupled cavities, the extension to 
bidimensional configurations is certainly possible. We iden- 
tify in semiconducting photonic crystal microcavities (doped 
with multiple impurities) a foreseeable potential setup for our 
proposal BSfl. By requiring a small number of impurities per 
cavity, our proposal demonstrates the interest of a working- 
point which is intermediate with respect to the simulation of 
the Bose-Hubbard model (expected when each cavity is off- 
resonantly interacting with a large number of two-level sys- 
tem) and the fully microscopic polaritonic dynamics 

The physical model- In order to clearly elucidate the salient 
features of the mechanism under consideration, we address 
the simple case of a linear array of mutually coupled micro- 
cavities. For the sake of simplicity, we describe the cavities 
as M equally spaced and independent optical resonators, each 
occupying a site of the array. Each cavity contains N iden- 
tical two-level systems whose ground (excited) state we de- 
note as \g) (|e)). The specific nature of the two-level sys- 
tems is immaterial for the following duscussion. Therefore, 
for easiness of language, we refer to them as atoms until we 
address the physical setup we suggest. The cavities are suffi- 
ciently close to each other to allow for their mutual coupling 
through evanescent-photon hopping (see Fig. Q] for a sketch). 
The cavity-coupling drops exponentially with the distance so 
that we consider only nearest-neighbor interactions. 

Each atom interacts with the respective cavity via the 
electric-dipole coupling. As the volume of each cavity can 
be made very small (~ with A c the wavelength of the 
cavity field), the coupling within each cavity can be large 
The free Hamiltonian of the whole system is Hq = 
Y^jLi (w c a]aj + u eg J2k=i ®kj )' w h ere ^ c is the frequency of 
the cavity fields, u eg is the transition frequency of the atoms, 
Oj is the creation operator of mode j and is the z-Pauli 
matrix for the A:-th atom interacting with the j-th cavity of the 
array [we use units such that H = 1 throughout the paper] . In 
the resonant condition lo c ~ oj eg and in the interaction picture 
with respect to Hq, the coupling Hamiltonian reads 
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FIG. 1: Simplified schematics of the setup and representation of a 
possible descrete solitonic solution of the dynamics associated with 
the cavity fields (see the body of the paper). Each resonator ideally 
represents a defect microcavity in a photonic crystal, the correspond- 
ing ensemble of two-level atoms being embodied by impurities. 



Here, we have introduced the Dicke lowering operator of the 
atoms within a given cavity Sj = ^2k=l ®kj (with °kj = 



\g) k j(e\ and Sj = SJ ), the photon-hopping strength J and 
the coupling rate fl between an ensemble and its respective 
cavity. To fix the ideas, these parameters are taken to be 
real and homogeneous across the array so that Jj = J and 
f2j = fi, Vj, although the generalization to a pattern of cou- 
plings is quite straightforward. The first term (and its h.c.) 
in Ti. describes the tunneling-like coupling between the cav- 
ities while the second term accounts for the in situ resonant 
exchange of excitations between the N elements of an ensem- 
ble and the corresponding cavity field (energy non-conserving 
terms have been discarded by invoking the rotating wave ap- 
proximation, valid for 57 <C w c ). 

The extended character of the Hilbert space of the overall 
system makes the goal of tackling the dynamics described by 
ThL formidable. However, we can exploit the collective behav- 
ior of each ensemble by utilizing the Holstein-Primakoff (HP) 
transformation, which maps a physical collective spin into an 
effective bosonic degree of freedom 1 10]. Assuming N >• 1, 
we introduce the HP creation and annihilation operators bj 

and bj (j — 1 , . . , M) as 

St = VNbtAj, Sr = VNAjbj, S*| = SjS^ - N/2, (2) 



with [bj,bj] — Sji. Here, the operator A, = J 1 — Pjbj/N 
is necessary in order to guarantee that the Dicke operators Sj 
and S? satisfy the proper SU(2) algebra. In terms of effective 
HP bosons, Eq. (Q} takes the form 
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h.c. (3) 



The enhancement of the intracavity coupling due to the col- 
lective nature of the ensemble is clearly retrieved. The non- 
linear nature of the interaction between the two bosons enter- 
ing the dynamics is encompassed by the operators Aj, which 
depends on the number of excitations in the state of the ef- 
fective HP particle. The strength of any nonlinear effect re- 



sults from a trade-off between (bjbj) and N. For our dis- 
cussion, it is useful to expand Aj in power-series as Aj — 

Tdlo (i-2?K2'i!) 2 ^w^~- Althou g h an Y order of this expan- 
sion could be taken, for clarity of presentation we consider the 
situation of a "mesoscopic" number of atoms per ensemble, 
i.e. N is such that Aj~l- b]bj/2N + 0{N- 2 ), Vj. Phys- 
ically, this implies that the number of implanted two-level 
atoms per cavity has to be large enough for the HP transforma- 
tion to be valid but sufficiently small not to blur any nonlinear 
effect. By introducing the truncated expressions for Aj's into 
Hh p we clearly identify the form the non-linear mechanism 
has, up to the considered order of approximation 
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bjbj&j + h.c.) 

(4) 

The structure of this Hamiltonian is interesting. The term 
il^/N (bjCLj + bcij) is quadratic in the quadrature variables of 
the two modes and is obviously linear. On the other hand, 
the term VL{b^bjaj + b^b'jaj)/2^/N encompasses the nonlin- 
ear character of the interaction that, noticeably, is not of the 
cross-Kerr form [6]. It is worth stressing that, while in the 
usual approach to coupled-cavity problems, the hopping term 
in Eq. (0]i is diagonalized by means of a canonical transfor- 
mation introducing collective cavity modes fljj] , here we keep 
its structure which is instrumental to the following discussion. 
We notice that, by assuming the conditions for adiabatic elim- 
ination of the atomic degrees of freedom, each bj would ef- 
fectively become proportional to <ij and Eq. would reduce 
to the Bose-Hubbard model QJLJ] . 

Dynamical equations.- It is straightforward to find the 
Heisenberg equations of motion for the two bosonic species 
involved in the dynamics. These read 



2VN 



id t cij — — J(a 1+ i + a,-i) + fiv Nbj ■=b\b i , 
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id t b 3 = QVNaj - ;r ^=(2&tS :/ a j + b^a]) 



(5) 



The coupled nature and the discreteness of these nonlinear 
equations make their direct solution a formidable problem. In 
order to tackle such a challenge, we work in a way to recon- 
duct Eqs. <(5j into a manageable form, exploiting some plausi- 
ble physical assumptions. First, we would like to eliminate the 
difficulties related to the non-commutativity of the operators 
in Eqs. (O. Now, if the operating temperature of the system is 
low, we can neglect the fluctuations of the operators bj , cal- 
culated over the state the atoms are into and concentrate just 
on their mean values. Furthermore, we can take each cavity as 
prepared in a large-amplitude coherent states. This allows us 
to replace the operators appearing in the dynamical equations 
of motion with a corresponding complex scalar function 111 ill . 



We therefore operate the substitutions a; — > a; , a\ 



■j\ and 

Pi, b] — > /?*, where each complex quantity is a function 
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of time. Eqs. (0 now become 

id t ai = -J(a t+1 + a;_i) + QVN ( A 



T^' with depending on 
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(6) 
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As anticipated, these equations are easily generalizable to 
the case of arbitrary-order expansion in the nonlinear oper- 
ator Aj. Our approach is designed in order to allow for 
such a generalization. So far, we retained the discreteness 
of the equations of motion. For a clear understanding of 
our results, it is more intuitive to proceed towards the con- 
tinuous limit where the function 7 depends on a continuous 
position-variable x along the linear array of cavities (here 
7 = a, (3 and their h.c). This simply means that the inter- 
distance d between two consecutive atom-cavity systems is 
much smaller than the scale represented by the wavelength 
A of the wave-like excitation propagating across the array. 
That is, by taking k = 2w/X, it must be kd <C 2tt so that 
7j — * j(x,t). In this way, the photon-hopping contribu- 
tion to the first of Eqs. (|4]i can be modified considering that 
ctj+i + Oj-i — > 2a(x,t) + d 2 d xx a(x,t). This changes 
the equation for a(x,t) into idta(x,t) = — 2Ja(x,t) — 
Jd 2 d xx a(x,t) + nVN[(3(x,t) - i)| 2 /3(x, i)/27V], re- 
minding of a Schrodinger equation with a cubic non-linearity 
(SCN) ||4|], which is known to admit solitonic solutions. How- 
ever, the equation at hand is still quite different from an SCN 
as the nonlinear term couples the dynamics of the real bosons 
to that of the effective ones. For compactness, we drop any 
explicit dependence on the variables x and t. 

Rather than using brute force numerical methods for the 
solution of the problem (such as Newtonian relaxation tech- 
niques lfl2lp . we would like to get a clear picture of the phys- 
ical process. We therefore decide to use an approach that is 
able to capture the essential features of the physics behind 
this problem. In Ref. 11311 . the so-called multiple-scale tech- 
nique has been used in order to investigate gap solitons in non- 
linear periodic structures. The flexibility of this technique and 
its successful application in II 1311 suggest its adaptation to the 
present situation. The method is based on the expansion of 
both time- and space-derivatives in terms of mutually inde- 
pendent time- and length-scales v p = p?v (v = t, x) ac- 
cording to d v = J2 p >o V p d Vp - Similarly, 7 = J2 p >i ^ p l (p \ 
where we have introduced the small parameter [i. In its 
essence, the multiple-scale method utilizes a separation of 
spatial and temporal scales in analogous to other standard 
techniques used in quantum optics such as the slowly vary- 
ing envelope approximation 16[]. Following 111311 . we stop the 
expansion at order p = 3 in 7, which is enough, in the con- 
ditions of small nonlinearity at hand, to encompass the im- 
portant aspects of the system's dynamics. We now proceed 
to solve the problem corresponding to a scale of order p — 1 
and use it into the one at order p. By replacing the scale- 
expansion into the equations of motion and collecting terms 
corresponding to the same power of /i, we get the univer- 
sal structures (id to + 2 J + Jd 2 d XoX „)a^ - fly/N^ = 



r (p) , id t0 pw - nVNa^ - x fj wiiu x a (i 

higher-order derivatives of cS p ~ k ^ and /JvP-'O (k %y The 
study is performed stepwise: At order /1 we have T^ 1 ^ = so 
that and ft 1 ' satisfy linear equations in the slow variables 
xq and to which are combined together to give 

{d toto - 2iJ3 t0 - iJd 2 d t0X0X0 + fl 2 N)a^ = (7) 

with /?W = -ifl^/N Ja^dt . In order to solve Eq. © 
we take the plane-wave ansatz cr 1 ' — g e l { kx o~^to) w jth 
£ an envelope function depending on space and time scales 
faster than xq and to and uj the frequency of such a car- 
rier excitation. The associated solvability condition leads 
to the dispersion relation ui±(k) = — J(l — d 2 k 2 /2) ± 
\/J 2 (l - d 2 k 2 /2) 2 + Q 2 N, thus defining effective "optical" 
and "acoustic" branches (corresponding to uj+{k) and uj-(k) 
respectively) reminding us of a linear diatomic crystal. The 
difference with our case is that the two species in the crystal 
(i.e. true and effective bosons) share the same site-location. 
We remark that uo±{k) are the long-wavelength approxima- 
tions of the branches coming from the exact dispersion rela- 
tion for the discrete equations of motion © in the linear limit. 
This reads w ± (fc) = -Jcos(fcd) ± y/ 'J 2 cos 2 (kd) + n 2 N , 
which becomes co± (k) for kd <C 1 and exhibits a band-gap at 
the edge of the first Brillouin-zone [150. Here, in order to be 
clear, we concentrate on the optical branch. 

The previous level of solution leaves £ unknown. Its form 
is determined by going to higher order in p and imposing the 
appropriate solvability conditions (such as the nullity of any 
secular term, as required in order for the multiple-scale ap- 
proach to hold lfl3ll ). At order /1 2 this results in the equation 
(d tl + v g+ d xi )£ = 0, so that the envelope function must de- 
pend on the variable £ = X\ — v g+ ti with the group velocity 
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J 9+ (ui 2 + +fl 2 N) 

finally leads to the generalized SCN 



-Wp. The iteration of such an approach 
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0. 
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where c x = (Jd 2 uj\ + S1 2 Nv 2 + )/[uj + {uj\ + Cl 2 N)], c 2 = 
2Q 4 N + n 2 N)}, £ = (l//u)e, \ = £/m and t 2 = 
H 2 t, as implied by the multiple-scale definition. Eq. ([8]) can 
be exactly solved by means of well-known inverse-scattering 
methods (ISM) and is known to admit the solitonic solution 
e = T]y/2c 1 /c 2 sech{r)[x - 2 Cl at] - v } e ™(x-^i(<y 2 -v 2 )t-t4>o 



for c\c 2 > H1411 . which is the case here. We have thus 
retrieved a solitary-wave behavior in the linear array being 
studied which was the central aim of our study. Here, ry, a, v 
and $0 are integration constants to be determined from the 
boundary conditions attributed to each specific physical prob- 
lem. The inclusion of higher-order terms in the expansion 
of Aj 's (implying the use of lesser atoms per ensemble) re- 
ft*) 

suits in the modification of and, consequently, of ci j2 . 
Multi-soliton solutions can also be found through the ISM 
theory [14]. The transition from solitonic dynamics to the 
linear regime occurs for c\ c 2 . In the long-wavelength 
approximation and in the realistic situation of J <g; 57, this 



4 



corresponds to sf~N ^> Q/J, which imposes limitations on 
the number of atoms per cavity, in agreement with our anal- 
ysis above. As an important remark, we notice that the con- 
tinuous limit adopted here for clarity of exposition can be by- 
passed by approaching the discrete dynamical equations so 
as to find a solution valid within the whole Brillouin zone. 
This can be done with minimal changes to our approach 
and without affecting the conclusions of our study by taking 
7/ = E P M p 7 (p) (£, t, fa) with £, - v{ld - At), r = /it, 0, 
representing a phase difference between each system and A to 
be found via the solvability conditions lllrSll . 

Physical setup.- We now briefly discuss a potential candi- 
date for the observation of the predicted phenomena, which 
we identify in an array of coupled photonic crystal microcav- 
ities QSQ. A sketch is provided in Fig.|2] We stress that the lin- 
ear arrangement addressed here can be generalized with some 
efforts to a bidimensional configuration by replacing the cav- 
ity adjacency-matrix, which determines the hopping term in 
the first of Eq. (0, with an appropriate tensor. Each cavity is 
formed as a single defect (created, for instance, every other 
site) in the lattice-pattern of the crystal. The spacing between 
the cavities can be made as small as tenth of the wavelength 
of the radiation confined in the crystal (which needs to be a 
slab of thickness ~ 0.5d in order to confine radiation on the 
plane la]). The microcavities are doped with two-level im- 
purities (for instance, substitutional Si donor impurities in a 
GaAs photonic crystal), N being determined by properly ad- 
justing the bulk doping density, N ~ 10 being experimentally 
realistic [9]. As stated, the small mode-volume of each cavity 
allows for large atom-field coupling rates. A proper pattern 
of the array, together with improved quality of the resonators, 
allow for D,/Jr~j 10. The effects of atomic spontaneous emis- 
sion and relatively small cavity quality-factor in the formation 
of solitary waves is beyond the scopes of this work and is the 
focus of ongoing investigation. 

Remarks.- We have addressed the dynamics of an array of 
coupled cavities from the perspective of intrinsic nonlinear 
optics. The nonlinearity arises from the coupling of the cavi- 
ties with ensembles of two-level atoms. Such built-in nonlin- 
earity and the periodic photon tunneling compete so as to cre- 
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FIG. 2: Sketch of the physical setup, (a): A photonic crystal with de- 
fect microcavities (formed every other site). The microcavity spac- 
ing is d. (b): Each cavity contains N impurities and is coupled to its 
nearest-neighbors via photon-hopping mechanism. 



ate solitary waves. We have used a simple iterative method to 
demonstrate such a behavior in a system of coupled photonic 
crystal microcavities. Our study contributes to the affirmation 
of systems of coupled cavities as flexible quantum simulators 
and represents an original approach to their dynamics, inter- 
estingly counterposed to current mainstream 
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